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We examine annihilation of light bosonic Dark Matter into pair of photons in model-independent 
way. We consider the simplest generic Lagrangian describing such process and then compare results 
to the available experimental data. Then we match our results with particular Dark matter models 
and determine possible constrains onto parameter space of those models. 



The presence of cold Dark Matter (DM) in the Uni- 
verse provides explanations to several observational 
puzzles and is an established fact nowadays. How- 
ever despite numerous experimental efforts the nature 
of DM remains a mystery. Many elementary parti- 
cle theories beyond the Standard Model (SM) have 
in their content at least one electrically neutral, sta- 
ble, weakly interacting particle. In order to make se- 
lection between those models a consideration of DM 
properties using different observables with the mini- 
mal number of assumptions is needed. In this paper 
we provide a model-independent approach for anni- 
hilation of light cold scalar dark matter. Such limi- 
tations are motivated by WMAP observations which 
ruled out warm DM [l| and the fact that Lee- Weinberg 
limit that forbids light dark matter can be avoided 
for non-fermionic Dark Matter particles Q. It means 
that independent constraints of such models are use- 
ful tool for discrimination of different theories. It 
can be argued that combined constraints from heavy 
quarkonium decays, astrophysical observation, and di- 
rect DM detection experiments can limit parameter 
space of such Dark Matter candidates Q. Here we 
report on astrophysical observations. 

Due to large available amount of data we used 
gamma ray flux as an experimental observable to con- 
strain properties of DM . We compare data from 
EGRET with theoretical calculations of flux from 
4>4> — * 77 process. This process is suppressed com- 
pared to the (jxj) —>■ Xj, however it provides very dis- 
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tinct spectrum feature and is very easy detectable. 
Current data from EGRET do not have any signs 
of monochromatic lines in photon spectrum, meaning 
that flux from (jxp — > 77 is below diffuse background. 
This fact can be used to derive constraints on proper- 
ties of Dark Matter. 

The paper is organized in the following way. We 
introduce generic model-independent lagrangian de- 
scribing DM-photon and compute gamma ray flux 
from DM annihilation. 



I. PHOTON FLUX FROM DM 
ANNIHILATION 

In general, annihilation process can be described 
with an effective Lagrangian of the following form: 

C eff = A X UF^F^ + AiMfF^F"" (1) 

Any additional operators will be of higher dimension 
and therefore their contribution will be suppressed. 
Completing the textbook level calculation of the cross- 
section one can get 



(\A^ + 2\A 2 \ 2 ) 



(2) 



Assuming that we deal with non-relativistic light 
dark matter, we present transferred energy as 



,2 ,., ..,2 ,2,., , f vi-^ n (3) 



s«4 Al 2 + (pl-p 2 ) 2 = (2/i) 2 (l + (- 



where fx is a mass of Dark Matter particle. Another 
assumption we make is that DM particles are 
distributed according to Maxwell-Boltzman 
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distribution. After we expand cross-section around 
s — 4/x 2 and average it with MB distribution 



S*/*(\A 1 \* + 2\A 2 \*) 



(av) 



A + B- 



6kT 



A + B{v{-V2)\4) 
(5) 



where A and B denote the following combinations of 
Wilson coefficients A 



A = p?(\A x \ 2 + 2\A 2 \ 2 ) 
B = A1 2 (^(|A| 2 + 2|A| 2 ) + 

r) A r)A 
f(2Re[A x ^]+me[A^])), 



(6) 
(7) 



ds 



with their values taken at s = 4/i 2 . It is worth 
pointing out that a fraction ~ ^ <C 1 which 
means that in most cases the contribution from the 
second term is negligible. 

The differential flux of photons produced by DM 
annihilations is 0] 



I(E^) 



dN~, 



{av 



dE 2p? 



(8) 



where -0 is the angle between the galactic center and 
the line of observation, and 



ds 



p 2 \r 



l.o.s 



An 



(9) 



is an integral along line of sight which depends on the 
choice of dark mater halo profile. As it was argued in 
6] , the maximum flux will be in the direction of the 
galactic center. The highest value of flux is given by 
the choice of Navarro- Frenk- White (NFW) profile for 
the DM distribution. This will provide us with the 
upper limit on the theoretical value of the photon 
flux and on the parameters of dark matter. Using 
results from the same paper we obtain: 



I(E,ip) = 7.3 x 10 _5< ^ 7 



dE 2p? 



-cm- 1 s- 1 sr- x GeV- 1 



(10) 

In this result the dependence on the particle physics 
dynamics is separated from the structure of Dark 
Matter halo. 

To proceed further we need to introduce mechanism 
of DM annihilation. As electrically neutral field DM 
can not be coupled to the photons directly. It is 
natural to assume that for the light Dark Matter 
(/i < 5GeV) the only relevant couplings are the ones 
that couple it to the Standard Model fermions. If we 
limit ourselves to the operators of the highest 
possible dimension six, the effective Lagrangian will 
take the form 



fWWV, 




(a) Direct channel 




(b) and a crossed one 
FIG. 1: Diagrams contributing to DM annihilation 



C = ^(CxOr 



c 2 o 2 ) 



(11) 



where A is a heavy mass scale, for example mass of 
heavy mediator that provides the interaction 
between the SM and the DM sectors. The operators 
are defined as 



0\ = mfc/xfri/jip 
2 = im,f(j)cj)tpj 5 ip 



(12) 
(13) 



The choice is such that they are hermitian and their 
Wilson coefficients Gj are real, if) are the SM fermion 
fields. There are only two types of diagrams that 
contribute to the annihilation process ( FigQ] ) 
For computation of annihilation rate the 
contribution from all possible fermions should be 
taken into account (i.e. summation over leptons and 
quarks performed). Summation is assumed and 
charge of loop fermions is denoted as Qf in the 
analytical expressions presented further in text. 
Explicit calculation of those diagrams using 
introduced generic lagrangian leads to the following 
Wilson coefficients A ;2 : 



16ClTO 2 7T 2 Q 2 



sA 2 



({Am) - s)C (0, 0, s, m}, m),m 2 f ) + 2) 



An 
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tC 2 rn 2 7r 2 (5 2 Co(0, 0, s, m 2 , m 2 , m 2 ) 



(14) 



/ 



Where, s is the Mandelstam variable 

s = (pi +P2) 2 = (ki+k 2 ) 2 

and 

Cob?, (Pi -P2) 2 ,P 2 ,TO 2 ,m 2 ,TO 2 ) = 
d 4 q 1 



itt 2 (q 2 - m()((q+ Pl ) 2 - mj)((q+p 2 ) 2 - mi) 
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is a Passarino-Veltman three-point function, and for 
particular set of parameters arising here it can be 
expressed analytically in the following form: 



C (0,0, 



s. m , m , m 



1 



tan _1 ( 



V47 



(15) 



It is worth pointing out, that our result for 
annihilation cross section via channel governed by 
operator Oi essentially reproduces the result of Q 
for annihilation of Higgs boson into two photons and 
does not vanish in the heavy fermion mass limit. 
Also, when the mass of the dark matter particle is 
close to the mass of the fermion in the loop, some low 
energy resonance states that increase annihilation 
cross-section might appear. Such a situation needs 
special treatment and is not considered here. 
Let us now use the experimental data to put 
constraints onto C\ and C%. Experimental data from 
EGRET can be parameterized in the following way 



I = I g ai+Ie X cm^s^sr^GeV- 1 



Igal 

where 

N (l,b) 

N (l,b) 



(7.32 ±0.34) x 10" 



E 



N (l,b) x 10~ 6 I -JL 1 



0.451GeV 

2.7 



-2.10±0.03 

(16) 



\GeVj 



y/1 + (1/35) Vl + [&/(!■! + 0.022|Z| 
for Ul > 30 



Since the highest flux will be from the direction of 
the galactic center, we need to compute flux at 
I, b = 0. There was no monochromatic peak observed 
at EGRET, which means that intensity of flux from 
dark matter annihilation is less than diffuse 
background. This leads to the following constraining 
condition: 



Ithe 



Upper Bound(I ex + I ga i) 



< 1 



(18) 



Assuming annihilation of DM particles that are at 
rest, photon spectrum will be a monochromatic line 
with _E 7 = fi. Detector measuring spectrum has 
finite resolution, thus instead of <5-function a 
spectrum integrated over some region of energies will 
be measured. Considering several different masses of 
DM particles leads to the following constrains that 




(a) Static DM 




002 0.004 0.OB6 A 2 



(b) Relative motion taken into account 

FIG. 2: Constrains on DM parameter space in case of 
static (a) and moving (b) DM. Filled regions are allowed 
parameter space for DM particles of different mass: red 
-fl — O.lGeV, green -/i = 0.5GeV, blue -/j, = lGeV, pink 
-/j, = 2GeV and yellow - fi= 5GeV 



can be placed on coupling constants 



-180° < I < 180° and - 90° < b < 90° 

0.5+ - ^ (17) 

Vl + (Z/35)Vl + (Vl-8) 2 
for \l\ < 30 

0.5 85 ' 5 



2.(m7 |'g 



C2 

lY 2 



14.02 -4 < 1 



for fj, = O.lGeV 

2 

+ 424.1 



for n = 0.5GeV 

2 



9i 

A 2 



v# < 1 



479.0 



477.6 



Ci 
A 2 

Ch 
A 2 



2396 



2057 



9i 

A 2 
A 2 



< 1 for fi = l.OGeV 



< 1 



for fj, = 2.0GeV 

2 

+ 64.55 



for fi = 5.0GeV 



A 2 



^ <1 



which are presented graphically on the Fig 2(a) 
However, in real-life situation the spectrum will be 
smeared due to thermal motion of DM particle, 
orbital motion of Earth, etc. and therefore energy 
spectrum will have a shape of a peak of finite width 
and finite height instead of (S-function. Assuming 
good enough resolution of experimental set up, this 
peak might be detected. Thermal velocity of DM 
particles is taken to be vdm = 9km/ s Q and taking 
into account orbital motion of Earth we approximate 
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spectrum as a Gaussian distribution with a = fiv/c 
where v = vom + v orbital- Such a choice of v^m 
among all experimental data provides us with the 
most narrow and high peak. It will be very easy 
detectable and will lead to the highest values of 
upper bounds on DM model parameters. The results 
after relative motion is considered are presented at 



EqfUand Fig 2(b) 



C 2 = (20) 
A = M h with M h > H5GeV 



Inserting these parameters into model-independent 
bounds derived in EqQ2j] and Eq[T5] leads to the 
following constrains onto parameters of this model: 



2 Ax • 10 1 I § 



4.12 x 10 5 



1.47 x 10 b 



,A 2 

for fi = 0.1GeV 
A 2 

for Li = 0.5GeV 
'Cx 
A 2 

for fi = l.OGeV 
1.04xl0 6 (g 
for [i = 2.0GeV 
2.61 x 10 4 f S 



1.36 x 10 5 f -p ) ; I 



1.84 x 10 6 



7.35 x 10 6 



9i 

A 2 
A 2 



< 1 



< 1 



(19) 



4.46 ■ \\t ( ) ; I 



A 2 

for /i = 5.0GeV 



8.86 xlO 4 -4 <1 



C2 

A 2 



Realistically, after taking into account all possible 
effects, contsrains will be somewhere between ones 



provided in EafT9l Fig 2(a) and ones given in Eafl9l 
Fig|2(b)1 



II. MODEL OF SM SINGLET SCALAR DM 
AS AN EXAMPLE 

As an example we consider DM annihilation in 
framework of Minimal Scalar Dark Matter model 
(see for example [l0(). In this model DM interaction 
with Standard Model fields is mediated by exchange 
of a Higgs boson. The model is restricted based on 
relic abundance calculations, however due to its 
simplicity it is perfect for testing of our approach. 
The matching conditions for the Wilson coefficients 
have the following form: 

Ci = A/2 



|A| < 16227 (—} for ^ = 0.1GeV 

|A| < 2715.3 for fi = 0.5GeV 

|A| < 1208.5 (t^) for M = l-OGeV (21) 

|A| < 1210.4 (^^j for n = 2.0GeV 

|A| < 6064 (^) for/i = 5.0GeF 



for case of static DM and 



|A| < 164 [ — I for u = 0.1GeV 
115/ 

I A / 41.2 [ ^ ) for fi = 0.5GeV 
115 , 

A ; 21 S ( ] i\ )V j, = jj)(7 ( V (22! 

|A| < 26.0 ( ^ ) for n = 2.0GeV 

1 1 

|A| < 163.7 ( — I for [i = b.OGeV 
115 / 



if we take thermal motion of Dark Matter particles 
into consideration. 

As one can see, the obtained constraints are not very 
restrictive for this particular model. However, 
consideration of the models with enhanced couplings 
(for example two Higgs doublet model) provides 
more strict constraints onto the parameters of the 
model 0]. 
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